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Abstract 

In this paper we present a complete next-to-leading order QCD calculation of 
the Xj( 3 -Pj ; J = 0, 1, 2) hadronic decay width. We include the NLO colour-octet 
contribution, as defined in the Bodwin, Braaten and Lepage formalism. We extract 
an estimate of the colour-octet parameter #8 from the charmonium decay data. 
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1. Heavy quarkonium (HQ) systems are among the most interesting objects that 
nature gave us to explore perturbative Quantum Chromodynamics (QCD). The predic- 
tions of their production cross sections and decay rates were among the most important 
tests of the early-time QCD. Nowadays there is a large renewed interest in the physics 
of heavy quarkonium, due above all to the recent discovery of the surprisingly big dis- 
crepancies between data and theory in the high-p^ charmonium cross-section production 
at the Tevatron [![]. Bodwin, Braaten and Lepage (BBL) recently developed a new 
formalism based on non-relativistic QCD (NRQCD) || to implement in a systematic way 
both the relativistic and the QCD corrections to the naive Colour Singlet Model (CSM) 
(for a recent review, see for example @]). This framework has been successfully applied 
to solve the charmonium anomaly at the Tevatron [?]. Another crucial test of pertur- 
bative QCD in the heavy quarkonium systems is given by the P-wave hadronic decay 
rates. In the CSM the inclusive hadronic decay rate of P-wave HQ states shows a sin- 
gular infrared behaviour || which is a clear signal that such a process is sensitive to at 
least another non-perturbative parameter beyond the usual wave function. In particular, 
the infrared problem of xj { 3 Pj) decay arises from the T(xj — > qqg) subprocess. The 
amplitude associated to this process diverges when the final gluon becomes soft. This 
ambiguity spoils the traditional factorization picture even at leading order in a s in the 
decay of the xi state into light hadrons (LH). In the BBL theory there is the solution of 
the xj decay problem. In the CSM, the heavy quark pair that participates in the hard 
annihilation process is in a colour singlet state and has the same quantum numbers as the 
physical bound state: the non-perturbative transition changes neither the colour nor the 
spin-parity of the heavy-quark pair. In this picture the xj decay occurs through the non- 
perturbative transition xj QQ^Pj 1 ^} (the upper right label indicates the colour state), 
which is parametrized by the derivative of the wave function |/?'(0)| 2 , followed by the 
annihilation of the QQ[ 3 Pj^} heavy quark pair. Bodwin, Braaten and Lepage suggested 
that the HQ wave function contains a non- negligible component in which the heavy quark 
pair is in a QQ[ 3 S[ 8 ^} state. This component leads to the HQ decay through the process 
QQ[ 3 S[ ] — > qq. The final state created by the colour-octet contribution is degenerate 
with the colour-singlet one in the kinematical region of soft final gluon. The colour-octet 
long-distance matrix element absorbs the infrared sensitivity of the colour- singlet term 
yielding an IR-finite result |§. In the NRQCD framework it is therefore possible to give 
a theoretical prediction of xj hadron decays avoiding infrared inconsistencies. 

Hadronic xj annihilation then gives a very important phenomenological test of the 
role of the colour-octet mechanism, and of the BBL theory in general. In ref. M Bod- 
win, Braaten and Lepage performed a phenomenological analysis of Xcj decays using the 
LO results for both the colour-singlet and the colour-octet contributions. They justified 
the neglect of the known NLO colour-singlet corrections with the observation that NLO 
accuracy would require inclusion of the yet unknown NLO colour-octet coefficients. In a 
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following paper [ID] it was argued that it is justified to include in the analysis the available 
QCD NLO colour-singlet terms, because the octet contribution does not depend on J even 
at NLO. Therefore the higher-order colour-octet coefficient can be simply reabsorbed in a 
redefinition of the LO octet wave function without changing the values of a s and |i?'(0)| 2 
extracted from the fit to the experimental data. 

In this work we perform the NLO calculation of the colour-octet coefficient complet- 
ing the picture of the xj hadron decay at order 0(a s 3 v 5 ) in the sense of the BBL double 
expansion. The knowledge of the NLO octet correction allows us to extract the value of 
the parameter H 8 using the results of ref JTU] . 



An analogous calculation relative to the newly discovered l P\ state [O] has recently ap- 
peared |L2 . 



2. The xj state is represented as a Fock-space vector superposition of heavy quark 
pair states of different spin, angular momentum and colour, possibly accompanied by 
gluons 0: 

\ X j) = 0(1)\QQ[ 3 PP\) + 0(v)\QQ[ 3 s[%) + ■■■, (1) 

where v is the relative velocity between the bound quarks. The first term of eq. ([!]) 
represents the conventional colour-singlet configuration and the second one corresponds 
to a colour-octet heavy quark pair in the QQ[ 3 <Sf^] state accompanied by a gluon. The 
small relative bound quark velocity v splits the physics of heavy quarkonium into two well 
separated energy scales, allowing a formal factorization of the physical observables into 
perturbative short-distance kernels describing annihilation of the heavy quark pair and 
soft non-perturbative coefficients. In the BBL factorization framework, the low energy 
heavy quarkonium physics is described by the NRQCD Lagrangian which has a physical 
ultraviolet cutoff A . The short- distance annihilation effects are implemented including 
4-fermion interactions in the Lagrangian: 

$ ^4—fermions = 1 l®ni^)i (2) 

m n 

where m is the mass of the heavy quark. Both the NRQCD operators O n and the short- 
distance coefficients f n depend on A, but their product does not. The operators O n have 
well-defined scaling rules with velocity v and the coefficients f n have a QCD perturbative 
definition. Equation (H) can be actually read as a double a s and v expansion. For our 
study, the relevant operators O n are: 
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O 1 ( 3 P ) = \tf (-iB • <r) ^ (-^D ■ a) ^ (3) 
dfPi) = ^ * g x a) ■ 0t *g x a) V (4) 



0i( 3 P 2 ) = (-~BMJ 0t (-~5 ( V')J i> (5) 
8 ( 3 S' 1 ) = V>VT a . 0t aT > (6) 



The xj hadronic decay width can be written as: 



r( Xj — > LH) = 2 hatfPj) (Xj|0l(3 f j) ' Xj) + 2 Im /8 ( 3 ^) ^Wfi)l^> (7) 



The short- distance coefficients can be extracted by matching NRQCD and full QCD 
amplitudes 0. The NRQCD matrix elements can be determined phenomenologically or 
calculated on the lattice. 

Defining: 

= (xj\Oi( 3 Pj)\xj) H 8 (A) = (^MjllAM . ( 8) 

m 4 m 2 



we can rewrite the xj width as follows 



T( X J -> = rjfPj^ -> L#)#i + f 8 ( 3 5} 8) -> LH)H 8 = fi(J)Hi + f 8 H 8 (9) 



Velocity and mass scaling of the matrix elements of the relevant NRQCD operators are 
H\ ~ mv 5 , H$ ~ mv 5 . The QCD leading-order colour-singlet short- distance coefficients 
are of 0(a 2 s ) for xo an d Xi an d of 0(a 3 ) for xi states. On the other hand, the QCD 
lowest order colour-octet short- distance process QQ[ 3 S[ 8 ^] — > qq is of 0(a 2 ) while the 
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process — ► gg with both gluons on the mass shell is forbidden. Therefore a 
consistent perturbative picture of the xi decay at order a s 3 v 5 requires the calculation 
of the QCD NLO colour-octet contribution. For ease of reference, we collect here the 
expression for the xj decay widths including the NLO colour- singlet terms || and the 
LO colour-octet terms [HI 



T{x^LH) = ^a 2 H 1 \l + ^C ] 



log f + H 8 



T(x2^LH) = l7ia s 2 H 1 [l+ ^C 2 U +n 



(10) 



+ n Lt |a s 2 



i^ilogf +H 8 



where 



^)C F + n Lf (-if + | log 2 



°2= (ff -S7r 2 -21og2)C A -4C F + n Lf (-il + |l 0g 2 



11) 



01 - I 54 288 71 J + 81 



2x 



n Lf is the number of light quarks: n Lf = 3 for charmonium and n Lf = 4 for bottomonium 
states and a s = a s (m). Hi is related to the derivative of the wave function through the 
relation: 



9 [^(Q) ! 2 

2tt m 4 



#i 



1 + Oiv 2 



(12) 



We denote by £ a momentum scale that regularizes the soft divergence associated to the 
Xj ~^ QQ9 process. £ was usually related to the binding energy of quarkonium. Notice 
that neither the colour-octet term nor the coefficient of the divergent logarithm depends 
on the quarkonium spin J; this fact makes a universal renormalization of the parameter 
H$ possible. In fact, considering only the universal piece (Ur) of xj widths we get: 
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Ur 



TV 

3 

7T 

n Lf — a s 
3 



16 a s ( m A 

Zi 7V A 



- #8 (6) 
7T 

n Lf -a s 2 i7 8 (m). 



(13) 
(14) 



The A-dependence of the colour-singlet coefficient is consistent with that specified by the 
RGE for H R HI. 



3. We perform the calculation of the full NLO QCD colour-octet contribution to the xj 
decay widths term using the dimensional regularization scheme to regularize UV, IR and 
collinear divergences. We work in D= 4 — 2e dimensions. If we define: 



,(0) 



1 - e 



ira. 



'A-k^V r(l-e) 



3 - 2e V M 2 J T(2 - 2e) 



(15) 



where M = 2 m, then the D-dim Born colour-octet short- distance coefficient assumes the 
form: 



^(Born) 



n f (0) 



(16) 



The NLO correction to T 8 consists of real and virtual emission of gluons. 



Real emission 

The real correction to the short-distance colour-octet xj annihilation term is repre- 
sented by the two processes QQ[ 3 S{ 8) ] -> ggg and QQ[ 3 S{ 8) ] -> qqg . The calculation of 



the D=4 QQ[ 3 Si 8 ' ) ] — > ggg amplitude can be obtained via crossing from the results of ref 
0. This amplitude is completely IR and collinear finite because the two-gluon leading 
order amplitude vanishes. The calculation of the three-gluon real contribution is therefore 
straightforward. We obtain []: 



ri 0) -5 ( — • — ~ 



a. 



TV 



73 67 



4 



36 



(17) 



lr rhe expression reported here assumes implicitly N c 
matrix elements is not reported in the result of ref. 



3, since the explicit N c dependence of the 
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On the contrary the process QQ[ 3 s[ 8) ] — ► qqg shows IR and collinear poles that we ex- 
pect will cancel when adding the virtual correction. The D-dimension QQ[ 3 Si 8 ^] — > qqg 
amplitude that we obtain is in agreement with ref in the D=4 limit. It leads to the 
following width contribution 



dT 



{qqg) 



647ra 



3™ 3 



Af 2 (3-2e) 



t 4- u + 2M 2 s - e(t + it)' 



tu (s - M 2 f 



d(PS)[qqg}(lS) 



where s = (q + q) 2 , t — (q + g) 2 , u = (q + g) 2 . Performing the Mandelstam variable sub- 
stitution s = M 2 (l — x), t — M 2 xy, u = M 2 x(l — y), the phase space assumes the follow- 
ing form: 



d{PS)[qqg] 



M 2 (An^ 2 \ 2e 



1 



i28^ 3 \ m 2 rri - e 



-x 



x 2 (l — x)y(l — y) dxdy. 



(19) 



Integrating over the phase space we get 



i(qqg) 



f (o) a » 

J- s 



riLfi's"'— /(e) 

7T 



1 3\ _ 1 _ /19 
+ — I + c A — + C F 



e 2 2e 



l 2e 



7T 



4 3 



11 



~-C A 



(20) 



where 



/(e) 



4:7TfI 



2\ 6 



r(i + c) 



(21) 



Virtual emission 

The diagrams contributing to NLO virtual emission are shown in fig. (1). In table 1 
we list the contribution T>k of each diagram with the relative colour factors. The virtual 
colour-octet width can be written as 



^(Virtual) 
8,QCD 



rf/( e p£,m 



n. 



(22) 



Summing all the virtual diagrams, we find 
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^(Virtual) 
S,QCD 



n 



f(0) 
Lf 1 8 



i + -m 

7T 



2bn 



+ -{c F -\c A 

v V 2 



A 



+ 



2e T 



C, 



l 2e T 



(23) 



where 



b = 
A = 

v = 



1 ( ll r 2 t 

— C4 ^Lf J- F 

2 V 6 3 



C F - 



1 - 



3 
4m 2 



-7T 



50 2, n 7T 2 
h - log 2 

9 3 & 4 
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n T , f T P 



(24) 
(25) 

(26) 



The Coulomb singularity disappears by performing the matching between the NLO full 
QCD and NRQCD amplitudes, yielding a finite result in the v — > limit. Summing the 
real and the virtual emission corrections, we obtain the a s NLO colour-octet decay width 
for xj states; we give the result for Cp =4/3, Ca = 3 and T F =1/2: 



^(NLO) 



^(Born) 



a 



■ MS {v) ( 107 3 



71 

a ,, , £t 

- TT^Lf + 46 log - — 
9 2 m 



-7^ + 21og2 + 

4 



+ r 



(0) 



a 



MS 



(/*) 



7T 



73 67 , 

1 7T* 

4 36 



(27) 



Choosing the renormalization scale /jl = m ( the same of the colour-singlet terms in 
eqs. (|i~0D ) we obtain the following expression of the a s NLO imaginary part of / 8 ( 3 S'i): 



Im / 8 ( 3 5!; 



6 V 



a K MS (m, 



n Lf +rz L 



a s MS (m) 



7T 



107 3 2 

— - 7 7r 2 -91og2 + 
6 4 



5 2 

-n Lf + -n Lf log 2 



a 



MS, 



ym) „ ( 73 67 9 

V 7 5 1 7T 2 

7T V 4 36 



(28) 
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Always keeping \i = m, we report below the numerical colour-octet corrections for char- 
monium and bottomonium: 



f (7VLO) a MS (m ) 

= 1 + 3.9 ^-pl [charm] (29) 

1 8 

i^ = 1 + 3 . 8 2^-p [bottom] (30) 



Using the results of ref. it is straightforward to obtain the best fit of the pa- 



rameter Hg for charmonium Xc including the NLO QCD effects, and we get the value 
H 8 {c) (m c ) = 3.1 ± 0.5 MeV. Taking m c = 1.5 GeV we obtain (x CJ |C 8 ( 3 ^i; rn c )\x c j) = 
(6.8 ± 1.1) x 10~ 3 GeV 3 . For completeness we recall that the fits of the other parameters 
obtained in |IJ are E x = 13.7 ± 2.3 MeV and a s m {m c ) = 0.286 ± 0.031. As discussed 
in the introduction these results are not affected by the inclusion of the NLO colour-octet 
corrections. Using the NRQCD scaling rules, we can obtain an estimate of the bottom 
octet matrix element H^ b \rrib) — 0.66 MeV. 

We now want to analyse the renormalization scale dependence of the NLO colour- 
octet decay widths compared with the leading-order ones. The results are shown in figs. 
2 and 3 for charmonium and bottomonium states, respectively. The normalization of 
the bottomonium width is achieved by using the estimate of the colour-octet parameter 
Hg- b \rrib) obtained above through the NRQCD scaling rules. For the running of two-loop 
a s we use the input A*^£ 5 = 160 MeV extracted from the fitted value of a s MS {m c ). 
The pictures show that the inclusion of NLO corrections significantly reduces the scale 
dependence of the processes. 

To conclude we notice that the calculation presented here can be used to compute the 
strong NLO qq — > QQf 3 ^ 8 ''] contribution to the total ip hadronic production cross section 
and the NLO colour-octet fragmentation function of the gluon into ip. Work on these 
issues is in progress. 
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Diag. V k 



fk 



2e uv ^ 26 
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2C F - \C A 



Table 1: Partial virtual QCD corrections to the process QQ[ 3 S[ ] — > qq 
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Figure 1: Virtual Feynman diagrams contributing to the process QQ^S^ } — > qq 
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Figure 2: Renormalization scale dependence of the colour-octet contribution to 
Xc hadronic decay width Tg = T$Hg 
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Figure 3: Same as fig. 2 but for Xb 
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